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Abstract 



A study of exclusive production of the DD system through initial-state radiation is performed in 
a search for charmonium states, where D = D° or D + . The D° mesons are reconstructed in the 
D — > K~ir + , D° — > K _ 7T + 7r°, and D° — > K~7r + 7r + 7r~ decay modes. The D + is reconstructed 
through the D + — > If ~7r + 7r + decay mode. The analysis makes use of an integrated luminosity of 
288.5 ftr 1 collected by the BABAR experiment. The DD mass spectrum shows a clear ^(3770) 
signal. Further structures appear in the 3.9 and 4.1 GeV/c 2 regions. No evidence is found for 
Y(4260) decays to DD, implying an upper limit g ^ (^eso)^ i- ) ^ ^ (95% confidence level). 
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1 Introduction 



The y(4260) structure was discovered in the initial-state radiation (ISR) process e + e~ — > 7/5_rvt + 7t~ J/ip 
at the energy where the total e + e~ hadronic cross section shows a local minimum [2]. Its spin- 
parity assignment, J PC = 1 is inferred because it can be produced by a single-photon annhilation 
mechanism. The measurement T(Y(4260) -> e+e~) x £(Y(4260) -> vr+vr- J/ip) = (5.5 ± l.Oi^) 
eV/c 2 and an upper limit estimate of the Y(4260) contribution to the total cross section gives a 
value for r, r7r j / /^(y(4260)) in excess of 1 MeV/c 2 [3], where T^ji^ denotes the partial width to 
the 7r + 7r~ J/ifi final state. In contrast, T^ w j/^ for the is (80 ± 33 ± 23) keV/c 2 [3], prompt- 

ing suggestions that the V(4260) could be an exotic object [Sj. The Y(4260) — > 7r + 7r~ J/-0 decay 
channel indicates cc content in the y(4260), which is energetically allowed to decay to pairs of 
charmed mesons. In this paper, we report a search for V(4260) — > DD using events accompanied 
by initial-state radiation at the T{AS) energy. 



2 Data selection 

We use data from the BABAR detector [H] at the PEP-II asymmetric-energy e + e~ storage rings, 
located at the Stanford Linear Accelerator Center (SLAC). These data represent an integrated 
luminosity of 288.5 fb" 1 collected at y/a = 10.58 GeV, near the peak of the T(4S) resonance, and 
approximately 40 MeV below this energy. 

Charged-particle momenta are measured in a tracking system consisting of a five-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer central drift chamber (DCH), both situated in a 
1.5-T axial magnetic field. An internally reflecting ring-imaging Cherenkov detector (DIRC) with 
fused-silica bar radiators provides charged-particle identification. Kaon and pion candidates are 
selected based on a likelihood calculated from the specific ionization in the DCH and SVT, and the 
Cherenkov angle measured in the DIRC. A CsI(Tl) electromagnetic calorimeter (EMC) is used to 
detect and identify photons and electrons. Neutral pions are reconstructed from pairs of photons, 
each depositing at least 100 MeV in the EMC and having an invariant mass between 115 and 
155 MeV/c 2 . Photon pairs are kinematically constrained to have the tt° mass; pairs having a fit 
probability greater than 1% are retained as 7r° candidates. 

Candidate e + e~ — > jisrDD events, where ^isr denotes an ISR photon, are reconstructed using 
four decay channels: 

e + e" D°D°, D° -> i^7r + , D° -> i^+vr" (1) 

e + e - ^ 1ISR D 0d°, D° -> ET-tt+tt , D° -► K+n' (2) 
e+e" -> <y IS R D°D°, D° -> ^"vr+vr+vr", D° -> i^+vr" (3) 
e+e" -> 77 5R D + D-,D + D~ -> K^^i:-. (4) 

Charge conjugate decays are implicitly included. 

All charged daughter tracks of D meson candidates must be well-reconstructed; however there 
may be one additional track having fewer than twelve DCH hits, having a transverse momentum 
below 0.1 GeV/c, or having a distance of closest approach to the interaction point of greater than 
1 cm in the plane transverse to the beam or 10 cm along the beam direction. Events having an 
extra 7r° candidate are vetoed. The tracks of each D candidate are topologically constrained to 
come from a common vertex; D candidates with a vertex fit probability greater than 0.1% are 
retained. Subsequently, each DD pair is refit to a common vertex with a beam spot constraint; 
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DD candidates with a fit probability of greater than 0.1% are retained. When calculating kinematic 
properties of DD pairs, the energy of the D mesons with daughter tracks that are all charged is 
recomputed using the nominal D mass value [7j. For candidate D° — ► K~tt + tt° decays, a one- 
constraint fit to the nominal D° mass is performed. Events having a DD pair with an invariant 
mass below 6.0 GeV/c 2 are retained. 

A distinguishing characteristic of e + e~ — ► jisrDD events is that the squared invariant mass 
of the recoil to the DD system (MM 2 ) is that of the initial-state photon. The peak centered 
on MM 2 = in Fig. summed over all four reconstructed channels, provides clear evidence for 
exclusive ISR production. The shaded histogram shows a background estimate derived from the 
two-dimensional D and D mass sidebands. 




MM 2 (D D) (GeV 2 /c 4 ) 

Figure 1: The recoil mass squared for ISR event candidates (MM 2 ) summed over all the four 
reconstructed final states. The shaded histogram is the background distribution estimated from D 
and D mass sidebands. 

Candidate events in the ISR region, defined as |MM 2 | < 1.0 GeV 2 /c 4 are retained. Observation 
of the ISR photon is not required since these photons are emitted preferentially along the beam 
directions; approximately 10% are reconstructed in the EMC. In Fig. |2J we show the mass of the D 
vs. that of the D, for each reconstructed channel and projections for each reconstructed D decay 
mode. 

The number of observed events for each reconstructed channel is given in Tabled The same 
table also shows the resolution from a single- Gaussian fit to the D-candidate mass spectra and 
the experimental DD mass resolution at the mass of the K(4260). The mass resolution for each 
channel agrees with values found for simulated e + e _ — * "fisRDD events. A ±2. 5a signal region is 
defined for each D channel, with sideband regions defined from -6.0 a to -3.5 a and 3.5 a to 6.0 a. 
Background estimates for each channel are inferred from the two-dimensional D and D sidebands. 
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Figure 2: Scatter plots for each channel and projections of the candidate D mass. 
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Table 1: Mass resolutions and event yields for each reconstructed channel. 



Channel 






D Mass Res. 
( MeV/c 2 ) 


Candidates 


Background 


Signal 


DD Mass Res. 
(MeV/c 2 ) 






\W -> K+TT- 


6.8 ± 0.2 


68 


4 


64 ± 8 


5.1 ± 0.1 






ir°,D° -> K + ir- 


13.0 ± 0.4 


121 


22 


99 ± 12 


5.2 ± 0.2 


D°^K~ 






5.0 ± 0.04 


111 


16 


95 ± 11 


4.6 ± 0.1 


D + -> K~ 


~7T~ 




5.6 ± 0.1 


116 


14 


102 ± 11 


4.6 ± 0.1 



BaBar 

preliminary 
200 




0.1 0.2 0.3 0.4 

Am(D° y) (GeV/c 2 ) 



0.1 0.2 0.3 0.4 
Am(D°7) (GeV/c 2 ) 



Figure 3: Spectra of Am(D j) = m(D 0, y) — m(D°) for (a) all events inside and (b) D° — > 
K~tt + ,D — > iT + 7r~ events outside the exclusive ISR region. The arrow indicates the expected 
position of D*°. 



A potential source of contamination in this sample of e + e~ — > jisrDD candidates could arise 
from the low MM 2 tail from e + e~ — ► jisrD^ D^*\ This background is partially suppressed by 
the extra-7r° veto; however radiative D*° decays or D*° — > tt°D decays with an unreconstructed 
photon daughter of the 7r° could survive. A direct search for D*° is made by plotting Am(D°j) = 
m(D°j) — m(D°) in Fig. Efa) for all D° and D° candidates within the ISR MM 2 region. For 
comparison, the same quantity is shown for D — > AT~7r + , D° — ► lf + 7r - candidates outside the ISR 
MM 2 region in Fig. Efb). A clear peak in Fig. Efb) corresponding to D*° — > is visible, but 
there is no evidence for any D*° contamination in the ISR sample. 
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3 DD mass spectrum 

The DD invariant mass spectrum for the ISR sample is shown by the data points in Fig. ^ A 
clear signal is seen for the -0(3770), which decays predominantly to DD. Additional enhancements 
coincide with the ^(4040), ■0(4160) and -0(4415) masses, which are observed in measurements of 
the total e + e~ hadronic cross section. Finally, a broad enhancement is evident near 3.9 GeV/c 2 . 
This structure has not been seen in the hadronic cross section measurements but is in qualitative 
agreement with coupled-channel model predictions [S]. 



60 




3.8 4 4.2 4.4 4.6 4.8 
m(D D) (GeV/c 2 ) 



Figure 4: The DD invariant mass spectrum, summed over all four reconstructed final states, with 
a fit that includes the y(4260) contribution. The arrow indicates the expected position of the 
F(4260). 



An unbinned maximum likelihood fit is performed using a signal shape described by four rel- 
ativistic P-wave Breit-Wigner distributions convoluted with a P-wave phase space function. The 
Breit-Wigner parameters are fixed to the values from a fit to the hadronic cross section |J2|. A 
Gaussian term is used to parameterize the enhancement near 3.9 GeV/c 2 . All Breit-Wigner and 
Gaussian terms are allowed to interfere by assigning them a free phase. A constant background is 
fixed to a value of 0.84 events per 20 MeV/c 2 bin as determined from a fit to the D and D mass 
sidebands. The fit yields a Gaussian-term mass of (3.909 ± 0.021) GeV/c 2 with a = (0.050 ± 0.007) 
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GeV/c 2 . The DD mass resolution at the Y(4260) mass, determined from Monte Carlo studies and 
shown in Table ^ is small compared to the widths of the fit structures and is neglected. 

We consider a possible Y(4260) contribution by performing a second fit, adding a fifth Breit- 
Wigner term with parameters determined from the Y(4260) observation in the J/i(m + Tr~ mode: 
m(Y (4260)) = 4.260 GeV/c 2 and r(Y(4260)) = 0.088 GeV/c 2 [I]. The Gaussian term is fixed to 
parameters from the first fit. The result, shown in Fig. EI yields a Y(4260) component with 7 ± 13 
events. The systematic uncertainty is evaluated as ±8 events by varying the mass and width of the 
four ip resonances, the Y(4260), and the Gaussian enhancement near 3.9 GeV/c 2 by one standard 
deviation, and by removing the phase space factor from the fit. All the resulting deviations are 
added in quadrature. 

The DD reconstruction efficiency for each channel, determined from Monte Carlo studies, in- 
creases with DD mass as shown in Fig. 
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Figure 5: Reconstruction efficiency for the four channels as defined in Eqs. (l)-(4) as a function of 
the DD invariant mass. 



Combining the fit result with world average D branching fractions, reconstruction efficiency, 
and the result from the BABAR , jisr^P 7T+7T ~ analysis yields 

£(F(4260) —* DD) _ 14±31 _ 



fi(Y(4260) -> J/V'vr+vr-) 



where statistical and systematic uncertainties have been added in quadrature. Because a statisti- 
cally significant signal for Y(4260) — > DD has not been established, this result can be recast as an 
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upper limit: 

■ < 7.6 at 95% C.L. (6) 
fi(Y(4260) -► J/^vr+vr-) v 7 

This limit is over an order of magnitude smaller than the value found for the -0(3770), another 
indication that the Y(4260) is not a conventional vector charmonium state. 



4 Conclusions 

In summary, ISR events have been used to explore DD production in e + e~ annhilation from charm 
threshold to 6 GeV/c 2 . We find evidence for DD decays of the V(3770), V(4040), 0(4160), and 
-0(4415), and a broad enhancement of DD production near 3.9 GeV/c 2 in qualitative agreement 
with coupled-channel model predictions. No statistically significant signal for F(4260) — ► DD is 
seen, leading to a 95% confidence level upper limit, g(y(42^^jj^n+n-) ^ 
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